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Abstract: We demonstrate dispersion tailoring by coupling the even and 
the odd modes in a line-defect photonic crystal waveguide. Coupling is 
determined ab-initio using group theory analysis, rather than by trial and 
error optimisation of the design parameters. A family of dispersion curves 
is generated by controlling a single geometrical parameter This concept is 
demonstrated experimentally on a 1.5mm-long waveguide with very good 
agreement with theory. 
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The linear and nonlinear propagation of waves depends crucially on the dispersion. The ability 
to control the dispersion is therefore a highly desired property of the physical medium sup- 
porting the propagation. For instance, the minimization of chromatic distortion in optical fibers 
im is the main concern in optical communications. Moreover, some nonlinear interactions e.g. 
parametric conversion/generation, involving different frequencies, are completely determined 
by the dispersion, through the phase matching condition||2|. 

In optical fibers and waveguides, dispersion control is possible because the waveguide disper- 
sion, which can be controlled through the geometry, can compensate for the material dispersion, 
which is material-dependent ||3]. In an optical waveguide, e.g. the so-called photonic wires, the 
dispersion can be controlled by changing the height and the width of the high index silicon 
wire on top of the silica layer |[3|. The ability to generate a prescribed dispersion function, also 
referred as dispersion engineering EUS), is a more demanding goal. By offering many degrees 
of freedom in the design, photonic crystals (PhC) fibers ||6| and waveguides Q [8] |9] |4] [TO) have 
proven to enable dispersion engineering. 

The particular dispersion of line-defect photonic crystal waveguides is often described as the 
result of an anti-crossing between index-guided (e.g. guided by total internal reflection) and 
gap-guided modes (which only exist in PhCs) 111 II . However the inherent complexity of pho- 
tonic crystals, as opposed to simpler optical structures, limits the understanding of the mecha- 
nisms underlying the control of the dispersion. As a result, in order to tune a PhC waveguide, the 
design is optimized regardless of the nature of the underlying mechanism, e.g. the changes in 
the mode field distribution or in the coupling between modes. This said, dispersion engineering 
has been demonstrated in PhC waveguides by optimizing the coupling of two waveguides||8]|2l 
or by coupling modes with different polarizations lfT2ll . 

In this paper we discuss a design procedure that, rather than optimizing some of the pa- 
rameters defining the PhC waveguide, defines a suitable perturbation which optimizes the anti- 
crossing of two, well-defined, modes. Especially breaking the natural symmetry of the waveg- 
uide in order to couple Odd and Even mode appears as an effective way to control the dispersion 
in PhCs |fT3][T4l[T5l ; still the symmetry should not be broken on a trial and error procedure. We 
demonstrate here that to obtain the desired dispersion features, depending on the initial design, 
the symmetry should be broken in a very specific way. This is illustrated by experimental re- 
sults. 

As an example, we consider here the even and odd modes, relative to the plane parallel to the 
propagation direction and perpendicular to the PhC slab. Those two modes are typically found 
in line-defect PhC waveguides and are uncoupled because they have different symmetries; how- 
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ever, the odd mode is rather regarded as a nuisance Umiting the spectral extent of single mode 
operation. Here, we show that, owing to even-odd mode coupling, a variety of useful dispersion 
profiles can be controlled by a single parameter, which is a quite desirable feature. 

To fix the ideas, let us consider first the case of two coupled PhC waveguides ||7l within the 
formalism of coupled mode theory (CMT). The dispersion of the super-mode resulting from the 
coupling Y{k) of the modes of the individual waveguides, with dispersion C0»,i,{i.2}{k)^ results 
from the diagonalisation of; 

Introducing the coupling results into a modification of the dispersion of the original modes, 
which can be controlled within some extent by playing with some parameters, such as the ra- 
dius of the holes or the separation between the two waveguides. We point out that y{k) could 
be modified, in principle, almost independently from fi)„,^,{,j / = {1,2}. 

Here, we explore the possibility of coupling the two modes of a single waveguide in order to 
control the dispersion. Operating with a single waveguide, rather than a set of coupled waveg- 
uides, has many practical advantages: the excitation (i.e. the coupling) of the mode is easier and 
the linear and nonlinear effective areas are smaller, resulting into stronger light-matter interac- 
tion, a desirable feature in nanophotonics. 

Specifically, we consider the odd mode, defined relative to the plane of symmetry x — z, where 
X is the direction of the propagation and z is perpendicular to the 2D PhC lattice. 

Let us consider a line-defect waveguide with dispersion represented in fig. [T];, calculated 
using the 3D plane-wave expansion method implemented in the MPB code|(16i|. The corre- 
sponding parameters, normalized with the lattice period, are: radius r/a = 0.26, slab thickness 
h/a ~ 0.41, width of the line defect W /a^ 0.95-^3, outward displacement of the first row of 
holes s/a~Q.\A, with radius reduced Xo r\/a — 0.23. The refractive index of the material is as- 
sumed to be 3.17. Here, the parameters have been chosen in order to lower the edge of the odd 
mode (at k = 0.5K, K being the reciprocal lattice vector); indeed, here the frequency spacing is 
only 8fa/c = 6x10^^. If a suitable coupling between these two modes, which are now closer 
in frequency, is introduced, then the anti-crossing will result into a change in the dispersion. As 
a consequence, controlling the coupling strength would enable, for example, the flattening of 
the lower frequency branch. 

Proceeding in the same way as for the coupled PhC waveguides, we now consider the cou- 
pling of even and the odd mode within the formalism of the coupled mode theory. The related 
eigenvalue problem is formulated by including the substrate mode which is mainly of the index- 
guided type. This is necessary because that mode might also interact with the other modes. The 
off-diagonal terms represent the coupling induced by a suitable modification of the PhC design. 
Namely the purpose here is to get a simple way to control the dispersion by tuning only the in- 
teraction between the odd and the even modes, without suffering from possible interactions 
with the other modes. 

(Oaddik) Joeik) Yos{k) 
'Yoe*(k) (Oeven{k) Yes{k) 
ts{k) fes{k) (>},ubstrate{k) 

Within the CMT formalism (which is a first order perturbation theory), the coupling between 
Jeo (k) is proportional to the spatial overlap between the odd and even modes, weighted by the 
perturbation Ae(r) of the dielectric function e(r) associated to the PhC structure, i.e.: 

Yoeik) = / / E:„,{r) ■ Ae(r) • dr (3) 
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In general, an arbitrary perturbation of the dielectric function e(r) would also introduce the 
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Fig. 1. PhC waveguide without (a) and witli (b) perturbation, c) Corresponding band dia- 
gram of the initial (perturbed) structure in black (red) line. The cyan markers indicate the 
modes for which the field distribution is shown d)-e). The thin blue line corresponds to the 
dispersion obtained for a different perturbation as discussed in Fig. [2] d-e) Spatial distribu- 
tion of the //--field, calculated at K=0.445, d) Even mode and Odd mode, e) even-like and 
odd-like mode. 



other diagonal and off-diagonal terms into eq. |2] 

Yooik) = J J EUiir) ■ Ae(r) dr 

Tee i^)= j j ^*even ' ' Eeven (?) dr 

Yssik) = j j El^bstrateir) ' Ae(r) • E„,hstrateir) dr 

%ld{r)-^£{r)-E,„harate{r) dr 
E*,.^,Xr) ■ Ae{r) ■ E,,,tsmite{r) dr 



(4) 



As a consequence, the perturbation would not only introduce even-odd mode coupling, but 
would also mix these modes with the substrate mode [yes, Yos)\ as well as with any other modes 
-which are not considered here-; and induce a frequency offset of the single dispersion profile 
through the diagonal terms. Hence the result of the perturbation is uncertain. If y^s = and 
Yos = 0, the dispersion is then easier to control, at least from a qualitative point of view. 

Hereafter we give a procedure to generate a perturbation Ae(r) such that Yeo{k) is non-zero, 
but also dominates also substantially over the other terms. Let us first look at the integrand 
E*^j{r) •A£(r) ■ Eeven{r)- The integral can be non-zero only if this is an even function. Con- 
sidering that the planar PhC waveguide belongs to the group of symmetry Cjv, that statement 
is equivalent to requiring the integrand to have at least the element of symmetry Ai (namely, 
the fully symmetric function). Conversely, if that condition is not respected, then integrand will 
have both positive and negative parts which would mutually cancel out. 

Proceeding similarly as in Ref. ifTTl . the Bloch modes in a PhC with a 2D hexagonal lattice 
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can be labeled according to the table[T| Based on the mode field distribution of the unperturbed 
structure in Fig. [TJl, we conclude that the even and the odd modes belong to the symmetry 
subgroups Bi and B2 respectively. The symmetry condition on the integrand in Yoe{k) is then 
formulated as; 

Ai^B2(E)X®Bi (5) 

where X is the symmetry of Ae(r), still to be determined. It is found that X must have the sym- 
metry A2. This implies that the dielectric perturbation is anti-symmetric with respect both to the 
x^O and the y = symmetry planes. This also implies that Yoo{k) = Yee{k) = 7v.s(fc) = Yse{k) = 
0, i.e., the only non-zero elements are Yeo{k) and Yso{k)- Indeed, since the even and the substrate 
modes have the same symmetry, they are not coupled together by a such perturbation. It can be 
observed that the coupling Yso{k) between the odd and the substrate modes is minimized if the 
initial frequency spacing between the two modes is large and if the changes of the dielectric 
function Ae(r) is localized where the field of the even and odd modes is much stronger than 
that of the substrate mode. This is in general the case in the region close to the line defect (see 

fig.IB). 

Based on that analysis, a modified PhC geometry is generated by translating the first row of 
holes by T /a = 0.1 along x in opposite directions, as is shown in fig.[TJ). It is straightforward 
to verify that the change introduced respects the A2 symmetry (fig.|2^). As a consequence, the 
original symmetry of the waveguide is now broken. The corresponding dispersion (figH];) now 
exhibits a kink near K = 0.445 and the lower (upper) band is translated to a lower (higher) 
frequency (spacing has now increased to dfa/c =1.5 10"^). The origin of the kink (i.e. zero 
group velocity dispersion points) is related to the change from a crossing behaviour (even and 
odd uncoupled) to an anti-crossing one, which implies that the modes repel each other 

We like to stress that there are several "odd" and "even" modes in PhC (to avoid confusion 
the modes should rather be labeled according the C2v group) which have different symmetry. 
As a consequence the way we broke the symmetry depends on the initial design. Hence in a 
different design fTOl, the dispersion is controlled by coupling the even mode with the substrate 
mode (which is initially closer in frequency); and this is done using a perturbation with the 
symmetry Z?i. In any case, the analysis based on group symmetry considerations gives access 
to the symmetry of the perturbation, and the location where it would be the most efficient. 
Therefore, we think this could help to reduce greatly the number of errors and trials done 
before eventually obtaining the desired dispersion. 

Because the waveguide symmetry is now broken, the strict classification of the modes as odd 
or even is not applicable. We therefore label the lower (upper) frequency mode as even-like 
{odd-like). The spatial distribution of the modes corresponding to these new bands is shown 
in fig. [T^ and clearly reveals the lack of a well defined parity (odd or even). The modification 
of the spatial distribution of the modes, i.e. their hybridisation through mode coupling, can be 
connected to the dispersion. 

In order to elucidate this point, we calculate the diagonal and off-diagonal terms introduced 




Fig. 2. a) Calculated \yee\, \yoo\ and \yeo\ coefficients for T=0.15a. Inset : corresponding 
Ae and its dominant symmetry (colors refers to opposite signs), b) Idem as a) but holes are 
moved in the same direction so the symmetry of Ae is now B\ instead of Ai- c) Group index 
dependence as a function of the wavelength calculated using MPB (resolution I8pts) for 
T= 0.05a, 0.1a, 0.2a. d) Corresponding measured dispersion in I.5mm-long waveguides. 
The lattice period is a = 465nm. 



in (IDi, which are shown in Fig. |2^. A closer look in the field distribution reveals that the per- 
turbation still keeps some features with the A\ symmetry and, as a result, small diagonal con- 
tributions and joo appear into the operator (|2]). We repeated the same calculation with a 
perturbation with equal magnitude but a different symmetry (namely B\), as shown in Fig. 
I2J3. Interestingly, while jee and jo,, are basically unchanged, the coupling between odd and even 
modes has disappeared. The dispersion (blue line on Fig.[T]: ) has no zero GVD points or equiv- 
alently, no anti-crossing is observed. This comparison confirms the minor role played by the 
diagonal contributions which are only responsible for a moderate steepening of the dispersion 
around K = 0.37, but do not create any inflexion point. 

It is now apparent that the features introduced in the dispersion are mainly related to the 
contribution above, which depends on T/a, a parameter defining the degree of asymmetry 
introduced in the structure. Therefore the change of the dispersion is controlled by an unique 
parameter. 

Let us consider now the dispersion, namely the group index as a function of the wavelength, 
as the parameter T is changed (fig.|2j;). Indeed, while a flat-band over more then 10«m is ob- 
tained with T=0.05, the other values generate a dispersion where the group dispersion changes 
sign twice and a local maximum of the group index increasing sharply as T is further increased. 

We have fabricated 1.5mm-long PhC waveguide on a 190nm thick air-suspended GalnP 
membrane (181 with the above specified designwhere T is varied between 0.05a and 0.20a. 
The use of mode adapters reduces the coupling loss below 1.5dB/facet [|19| . The dispersion 
is measured using an interferometric technique |20|. This technique also allow investigating 
the effects due to the disorder (multiple scattering etc.); these effects were found negligible in 
our waveguides. The result is shown is Fig. |2}l. The agreement with the predicted dispersion 
is very good. Residual differences (e.g. for T /a — 0.05) can be attributed to the fabrication 
tolerances (e.g. the radius of the holes). The important point is that the transition from one of 
these dispersion profiles to another follows continuously the change in the parameter T; which 
can be further adjusted in order to reach the desired dispersion profile. 

An important practical aspect in dispersion-engineered waveguides is their robustness 
against disorder, particularly their ability in the regime of low group velocity. For instance, 
a group index of 40 is reached for T /a = 0.2 (fig. |2};-d). Thus a group delay of about 200ps 
was obtained on a 1.5 mm long waveguide. This is quite a respectable target, comparable to 



the state of the art for delay Hnes based on nanophotonicsll2ll l22l l23l l24ll . Furthermore, in 
this slow light regime, the related increase of propagation losses is limited to about 5dB. For 
example, this allowed demonstrating efficient Four Wave Mixing in PhCs 125)1 . 

In conclusion, we have revisited the concept of dispersion tailoring based on mode coupling 
and we have implemented it in a single-mode line-defect PhC waveguide. The group theory 
analysis was used to determine a suitable perturbation of the PhC structure such that only the 
even and the odd modes are coupled. The change in the dispersion, and, particularly, the onset 
of two zero group dispersion points, is explained through the anti-crossing related to the even- 
odd mode coupling. A variety of dispersion profiles, including flat-band, is controlled through 
a single geometrical parameter, representing the degree of asymmetry which is introduced into 
the PhC structure. The same concept can also be extended to the coupling between two arbitrary 
modes (not only odd and even ones). This result is relevant to a broad range of applications, 
particularly parametric effects (e.g. Four- Wave Mixing) and it enables the optimization of the 
phase matching. We believe that it could be applied to other systems as well (e.g. PCF) and 
enable an even more accurate control of the dispersion. 

This work has been supported by the 7th framework of the European commission through 
the GOSPEL project (www.gospel-project.eu). We thank M. Santagiustina, G. Eisenstein and 
S. Trillo for enlightening discussions. 



